LOCAL STRONG SOLUTION TO THE COMPRESSIBLE 
VISCOELASTIC FLUID WITH LARGE DATA 



XIANPENG HU AND DEHUA WANG 

Abstract. The existence and uniqueness of local in time strong solution with large 
initial data for the three-dimensional compressible viscoelastic fluid is established. The 
strong solution has weaker regularity than the classical solution. The Lax-Milgram the- 
orem and the Schauder-Tychonoff fixed-point argument are applied. 



1. Introduction 

Elastic solids and viscous fluids are two extremes of material behavior. Viscoelastic 
fluids show intermediate behavior with some remarkable phenomena due to their "elastic" 
nature. These fluids exhibit a combination of both fluid and solid characteristics, keep 
memory of their past deformations, and their behaviour is a function of these old defor- 
mations. Viscoelastic fluids have a wide range of applications and hence have received a 
great deal of interest. Examples and applications of viscoelastic fluids include from oil, 
liquid polymers, mucus, liquid soap, toothpaste, clay, ceramics, gels, some types of suspen- 
sions, to bioactive fluids, coatings and drug delivery systems for controlled drug release, 
scaffolds for tissue engineering, and viscoelastic blood fluid flow past valves; see [6jlTTjl32"] 
for more applications. For the viscoelastic materials, the competition between the kinetic 
energy and the internal elastic energy through the special transport properties of their re- 
spective internal elastic variables makes the materials more untractable in understanding 
their behavior, since any distortion of microstructures, patterns or configurations in the 
dynamical flow will involve the deformation tensor. For classical simple fluids, the internal 
energy can be determined solely by the determinant of the deformation tensor; however, 
the internal energy of complex fluids carries all the information of the deformation tensor. 
The interaction between the microscopic elastic properties and the macroscopic fluid mo- 
tions leads to the rich and complicated rheological phenomena in viscoelastic fluids, and 
also causes formidable analytic and numerical challenges in mathematical analysis. The 
equations of the compressible viscoelastic fluids of Oldroyd type ([25, 26J) in three spatial 
dimensions take the following form [9| 120} [27]: 

Pt + div(H = 0, (1.1a) 
(pu) t + div (pu <g> u) - fiAu - (A + ^)Vdivu + VP(p) = div(/>FF T ), (1.1b) 
F t + u-VF = VuF, (1.1c) 
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where p stands for the density, u E R 3 the velocity, F E M 3x3 the deformation gradient, 
and P(p) the pressure which is a strictly increasing convex function of the density. The 
notation M 3x3 means the set of all 3 x 3 matrixes. The viscosity coefficients p and 
A satisfy the conditions that 3/i + 2A > and /x > 0, which ensure that the operator 
—pAu — (A + /x)Vdivu is a strongly elliptic operator. The symbol <8> denotes the Kronecker 
tensor product, F T means the transpose matrix of F, and the notation u- VF is understood 
to be (u • V)F. Usually, we refer the equation (jl.la[) as the continuity equation. 
We are interested in the Cauchy problem of (jl.ip with the initial condition: 

(p, u, F)(0, x) = (p , u , F )(x), x E R 3 - (1.2) 

The aim of this paper is to establish the local existence and uniqueness of strong solution to 
system (II. ip with large initial data in the three dimensional space R 3 . By a strong solution, 
we mean a triplet (p,u,F) with u(t, •) E W 2 ' q and (p(t, -),F(t, •)) E W 1 ' 9 , q > 3 satisfying 
(jl.ip almost everywhere with the initial condition (|1.2p . There have been many studies and 
rich results in the literature for the global existence of classical solutions (namely in H 3 or 
other functional spaces with much higher regularity) for the corresponding incompressible 
viscoelastic fluids, see [H El QUI US CEU HU QU EUJ E] and the references therein. For the 
compressible viscoelastic fluids (11.11) , the global existence of classical solutions in H 3 with 
small perturbation near its equilibrium for (II. lh without the pressure term was studied in 
|15j . and a local existence of strong solution near the equilibrium and a series of uniform 
estimates were obtained in [7] . One of the main difficulties in proving the global existence 
is the lacking of the dissipative estimates for the deformation gradient and the gradient of 
the density. To overcome this difficulty, for incompressible cases, authors in [14] introduced 
an auxiliary function to obtain the dissipative estimate for the classical solutions, while 
authors in [21] directly deal with the quantities such as Au + divF. Since we are concerned 
with the strong solutions in W 2,q which have weaker regularity than the classical solution 
in H 3 , we first linearize (jl.ip and use the Lax-Milgram theorem to obtain the solution 
to the linearized system, then we apply the Schauder-Tychonoff fixed point theorem to 
obtain the strong solution of (jl.ip . Using the standard notations W s,q (R 3 ) (H S (R 3 ) if 
q = 2) for the Sobolev spaces and setting Qt ■= [0, T] x M 3 for any T > 0, we can state 
our result on the existence and uniqueness as follows. 

Theorem 1.1. Assume that 

Po eW 1 ' q (R 3 )r\H 1 (R 3 ), u £H 2 (R 3 ), F €W 1 ' q (R 3 )nH 1 (R 3 ), 

for some q E (3, 6], and for some positive constants a, (3, and ro, 

en < po < P, ||uo||h 2 + IIPollw^rLH" 1 + ll F o||vyi.<jrLH" 1 < r o- 

Then there are positive constants T = T(ro), a(T, ro,a) ; and (3(T, ro,/3), such that, 
the Cauchy problem (jl.ip - (jl.2p has a unique strong solution (p,u,F) defined for (t,x) E 
(0,T) x R 3 , satisfying 

a(T,r ,a) < p< /3(T, r , /?); 
p E L°°(0,T; W 1,q (R 3 ) D H 1 ^ 3 )) n L°°(Qy); d t p E L°°(0, T; L q (R 3 )); 

u E L 2 (0,T; W 2 ' q (R 3 ) n H 2 (R 3 )); d t u E L 2 (0,T; H 1 (R 3 )); 
F E L°°{0,T;W^ q (R 3 ) DH^R 3 )) nL°°(Q ? ); d t F E L°°(0,T; L q (R 3 )). 
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The viscoelasticity system (jl.ip can be regarded as a combination between the com- 
pressible Navier-Stokes equation with the source term div(pFF T ) and the equation (|l.lc|) . 
As for the global existence of classical solutions of the small perturbation near an equilib- 
rium for compressible Navier-Stokes equations, we refer the interested reader to |22l [23] 
and the references cited therein. The global existence of strong solutions with small per- 
turbations near an equilibrium for compressible Navier-Stokes equations was also discussed 
in [231 [29]. Unlike the Navier-Stokes equations ([31 El IE]) , the global existence of weak 
solutions to (jl.ip with large initial data is still an outstanding open question. In this 
direction for the incompressible viscoelasticity, when the contribution of the strain rate 
(symmetric part of Vu) in the constitutive equation is neglected, Lions and Masmoudi in 
[19] proved the global existence of weak solutions with large initial data for the Oldroyd 
model. Also Lin, Liu, and Zhang showed in [12] the existence of global weak solutions 
with large initial data for the incompressible viscoelasticity if the velocity satisfies the 
Lipschitz condition. When dealing with the global existence of weak solutions with large 
data in the compressible case, among all of difficulties, the rapid oscillation of the density 
and the non-compatibility between the quadratic form and the weak convergence are of 
the main issues. For the inviscid elastodynamics, see [H [30l 13 lj and their references on 
the finite-time blow up and global existence of classical solutions. 

This paper will be organized as follows. In Section 2, we will recall briefly the com- 
pressible viscoelastic system from some fundamental mechanical theory. In Section 3 and 
Section 4, we will give the proof of the main theorem (Theorem II. ip . More precisely, Sec- 
tion 3 is devoted to the local existence of the system (jl.ip by the Lax-Milgram theorem 
and the Schauder-Tychonoff fixed-point argument, while the main goal of Section 4 is to 
prove the uniqueness of the solution obtained in Section 3. 



2. Mechanical Background of Viscoelasticity 

In this section, we discuss some background of viscoelasticity and some basic properties 
of system (jl.ip . First, we recall the definition of the deformation gradient F. The dynamics 
of any mechanical problem with a velocity field u(x,t) can be described by the flow map, 
a time dependent family of orientation preserving diffeomorphisms x(t,X), < t < T: 

f ix(t,X) = u(t,x(t,X)), 

} x(0,X) =X. [ ' ' 



The material point (labelling) X in the reference configuration is deformed to the spatial 
position x(t, X) at time t, which is the observer's coordinate. 

The deformation gradient F is used to describe the changing of any configuration, am- 
plification or pattern during the dynamical process, which is defined as 

. f)x 

F(t,X) = —(t,X). 

Notice that this quantity is defined in the Lagrangian material coordinate. Obviously it 
satisfies the following rule, by changing the order of the differentiation: 

9F(t,X) _ du(t,x(t,X)) 

dt - dx [Z - Z) 
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In the Eulerian coordinate, the corresponding deformation gradient F(t, x) will be defined 
as F(t,x(t,X)) = F(t,X). The equation together with the chain rule and (pTL|) . 

yield the following equation: 

d t F(t, x(t, X))+u- VF(t, x(t, X)) = d t F(t, x(t, X)) + gF(t '^' X)) • 

_ dF(t,X) _ du(t,x(t,X)) 
~ di ~ dX 

_ du(t,x(t,X)) dx 
~ dx dX 

du(t,x(t,X)) ~ u 
= F(t,X) = Vu-F, 

which is exactly the equation (|1 . lc|) . Here, and in what follows, we use the conventional 
notations: 

du ■ dF ■ ■ 

(Vu)ij = (Vu F)ij = (Vu) ik F kj , (u • VF)y = U k-Q^, 

and summation over repeated indices will always be well understood. In viscoelasticity, 
(jl.lcj) can also be interpreted as the consistency of the flow maps generated by the velocity 
field u and the deformation gradient F. 

The difference between fluids and solids lies in the fact that in fluids, such as Navier- 
Stokes equations [23], the internal energy can be determined solely by the determinant part 
of F (equivalently the density p, and hence, (jl.lcp can be disregarded) and in elasticity, 
the energy depends on all information of F. 

In the continuum physics, if we assume that the material is homogeneous, the conser- 
vation laws of mass and of momentum become [31 [T^l [3T] 

dtp + div(pu) = 0, (2.3) 

and 

d t (pu) + div(/m ® u) - pAu - (A + p)Vdivu + VP(p) = div((det F) -1 ^), (2.4) 
where 

pdetF = l, (2.5) 

and 

dW 

%(F) = Wj . (2.6) 

Here S, pSF T , W(F) denote Piola-Kirchhoff stress, Cauchy stress and the elastic energy 
of the material respectively. Recall that the condition (|2.6p implies that the material is 
called hyperelastic [20]. In the case of Hookean (linear) elasticity [T3 | IT4" 1 [16]. 

W(F) = \\?\ 2 = \tr(FF T ), (2.7) 

where the notation tr stands for the trace operator of a matrix, and hence, 

S{F) = F. (2.8) 
Combining equations (|2.ip - (|2.8p together, we obtain system (jl.ip . 



LOCAL STRONG SOLUTION TO COMPRESSIBLE VISCOELASTIC FLUID 5 

If the viscoelastic system (jl.ip satisfies div(poFo) = 0, it is verified in [21] (see Propo- 
sition 3.1) that this condition will insist in time, that is, 

div(p(t)F(t) T ) = 0, for t > 0. (2.9) 

Another hidden, but important, property of the viscoelastic fluids (jl.ip is concerned 
with the curl of the deformation gradient (for the incompressible case, see |13j I14j). Ac- 
tually, the following lemma says that the curl of the deformation gradient is of higher 
order. 

Lemma 2.1. Assume that (|l.lc|) is satisfied and (u, F) is the solution of the system (jl.ip . 
Then the following identity 

FlkVfoj = FijV^k (2.10) 
holds for all time t > if it initially satisfies (|2.10p . 

Again, here, the standard summation notation over the repeated index is adopted. 
Proof. First, we establish the evolution equation for the equality 

Indeed, by the equation (jl.lcp . we have 

dtViFij + u • VVjFy + V z u • VFy = V m u;V z F mi + V ; V 

Thus, 

F lk (dtViFij + u • VViFij) + F lk Viu ■ VF^- = Fi k V m UiViF mj + F lk ViV m \iiF mj . (2.11) 
Also, from (jl.lcp . we obtain 

ViF i:j (d t F lk + u • VF lk ) = V;F ij V m u i F mfc . (2.12) 
Now, adding OTTTT) and (j2TT2]) . we deduce that 

d t (Fi k ViFij) + u • V(F^V;Fy) = -F^V^u • VF^- + F ifc V m UjV;F mi 

+ F^V;V m UjF m j + ViFijV m UiF mk (2.13) 
m UjV;F m j + F;^V;V m UjF m j. 

Here, we used the identity which is derived by interchanging the roles of indices I and m: 
FzfcVjU • VFij = F Zfc V;U m V m Fy = ViFijV m uiF mk . 

Similarly, one has 

dtiFtjViFik) + u • V(F y V,F ifc ) = F y V mUi VjF mfe + FyVjV m UjF mii .. (2.14) 

Subtracting (l2~14l from ([2~13]) yields 

dt{FikViFij - FyVjFifc) + u • V(F, fc VjF fj - FyV,^) 

= V m Uj(F /fc ViF m j - F/jV z F mfc ) + V;V m Uj(F mj F/ fe - F m feFy). 

Due to the fact 

V z V m Uj = V m V z Uj 
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Hence, 



(2.16) 



in the sense of distributions, we have, again by interchanging the roles of indices I and m, 
V«V m Uj(F m jF/fc — F m kFij) = V;V m UjF m jF^ — V {V m \iiF mk Fij 

= V/V m UjF mj F^ — V m V;UjF^F m j 
= (V;V m Uj - V m V/Ui)F^F mj = 0. 
From this identity, equation (|2.15|) can be simplified as 

d t (FikViFij - FijViFik) + u • V(Fz fe VzFy - FyVjFifc) 
= V m Ui(F; fc V;F m j - FijViF mk ). 
Multiplying (l2~T6ll by F/ fc V/F„ - F^V/F^, we get 

d t \Fi k ViFij - FijViF ik \ 2 + u • V|FjjfcVjFy - F^-V/F^ 2 

= 2(F /fc V i F ij - - F /j V / F jfe )V m u,(F Zfc V i F mi - FyVjF^*) (2.17) 
< 2||Vu|| Loo(R3) X 2 , 
where M. is defined as 

M = max{|F/ fc V;Fj J - - FijViF ik \ 2 }. 

i,j,k 

Hence, (|2.17|) implies 

d t M + u • VM < 2||Vu|| i oo( R 3)A / i. (2.18) 
On the other hand, the characteristics of dtf + u • V/ = is given by 

-^X(s) = u(s,X(s)), X(t) = x. 

Hence, (|2.17p can be rewritten as 

^ < B(t, y)U, U(0, y) = M (y), (2.19) 

where 

U(t,y) = M(t,X(t,x)), B(t,y) = 2|| Vu|| Xoo(R3) (t, X(t, y)). 
The differential inequality (|2. 19[) implies that 

U(t,y)<U(0)exp(J B(s,y)ds 

M(t,x) < 7W(0)exp (J 2||Vu|| L oo (R 3 ) (s)^ . 

Hence, if A4(0) = 0, then A4{t) = for all t > 0, and the proof of the lemma is complete. 

□ 

Remark 2.1. Lemma 12.11 can be interpreted from the physical viewpoint as follows: for- 
mally, the fact that the Lagrangian derivatives commute and the definition of the defor- 
mation gradient imply 

~ d 2 Xj _ d 2 x,i _ ~ 

° Xk 13 ~ dx k dXj ~ dXjdx k ~ ° x > ik ' 

which is equivalent to, in the Eulerian coordinates, 

FifcViFy(t, x{t, X)) = FyVjFifcCt, x(t, X)), 
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that is, 

Fi k ViFij(t,x) = FijViF ik (t,x). 

Finally, if the density p is a constant, the equations of the incompressible viscoelastic 
fluids have the following form (see [2j [T2], [T3J, H31 USl [21] and references therein): 

divu = 0, 

<9iU + u- Vu-/iAu + VP = div(FF T ), (2.20) 
d t F + u- VF = VuF. 



3. Local Existence 

In this section, we will prove the existence part in Theorem I l.li The proof will proceed 
through four steps by combining the Lax-Milgram theorem and a fixed point argument. 
To this end, we consider first the linearized problem. 

Set 

W = {i> G (L 2 (0,T;H 2 (R 3 ))) 3 : G L 2 (Q T )} 
with the natural norm ||^>||w>, an d for q 6 (3,6] we define 
W = Wn (L 2 (0, T; W 2 ' q {R 3 )) n L°°(0, T; H 2 (R 3 ))f 

n {V G (L 2 (0,T;tf 2 (IR 3 ))) 3 : G L°°(0, T; L 2 (M 3 )), G L 2 (Q T ),^(0) = u }. 
Consider the following linearized problem: 

dtp + div(pi;) = 0, (3.1a) 
pd t u - pAu - (p + A)Vdivu = -pv ■ Vv — VP + div(/)FF T ), (3.1b) 
5 t F + f • VF = VvF. (3.1c) 
with the given v G W and the initial condition (jl.2p . 

3.1. Solvability of the density with a fixed velocity. Let Aj(x,t), j = l,...,n, be 
symmetric mxm matrices, B(x, t) an m x m matrix, /(x, t) and Vq(x) two m-dimensional 
vector functions defined in IR n x (0,T) and R n , respectively. 
For the Cauchy problem of the linear system on V G R m : 

n 

d t V + J2A j (x,t)8 Xj V + B(x,t)V = f(x,t), 

j=i I 3 - 2 ) 

V(x,0) = V (x), 
we have 

Lemma 3.1. Assume that 

Aj £ [C(0,T;F(R n ))nC 1 (0,T;ff s - 1 (R n ))] mxin , j = l,...,n, 

B G C7((0,r),F s - 1 (M n )) mxm , /€C((0,T),ff a (K n )) TO , V G tf s (M n ) m , 
wii/i s > 5 + 1 an integer. Then there exists a unique solution to (|3.2p . i.e, a function 

v g [C{[o,T),H s (R n ))nc 1 ((o,T),H s - 1 (R n ))] m 
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satisfying (|3.2p pointwise (i.e. in the classical sense). 

Proof. This lemma is a direct consequence of Theorem 2.16 in [23] with Ao(x,t) = I. □ 
To solve the density with respect to the velocity, we have 

Lemma 3.2. Under the same conditions as Theorem \l.l\ there is a unique strictly positive 
function 

p := S{v) E W 1>2 (0, T; L q (R 3 ) n L 2 (M 3 )) n L°°(0,T; W 1>q (SL 3 ) n H 1 ^)) 

which satisfy the continuity equation (|3.1ap . Moreover, the density satisfies the following 
estimate: 

ll V Plli<x.(o,T;M(K3))n^(K3) < (\\Vf*\\h(it*)nV<p*) + VT\\v\\w) exp (cVt\\v\\ w 

Here, and in what follows, the notation C stands for a generic positive constant, and 
in some cases, we will specify its dependence on parameters by the notation C(-); and 
W^ 2 (0,T;X) = {/:/,/«€ L 2 (0,T,X)}. In many cases below, we drop the constant C 
for the simplicity of notations. 

Proof. The proof of the first part of this lemma is similar to that of Lemma [3. 31 below, and 
can also be found in Theorem 9.3 in [24] . The positivity of the density follows directly 
from the observation: by writing (|3.1aj) along characteristics 4?X(t) = v, 

j t p(t, X(t)) = -p(t, X(t))dwv(t, X(t)), X(0) = x, 
and with the help of GronwalFs inequality, 

aexp(-Vf\\v\\ w ) < (infpo)exp (- J \\dxvv(t)\\ L co^ds\ < p(t,x) 

< (sup/) )exp (^J \\Aivv(t)\\ L oo^ds\ < /3exp(v / T||u|| w ). 

Now, we can assume that the continuity equation holds pointwise in the following form 

dtp + p divv + v ■ Vp = 0. 

Taking the gradient in both sides of the above identity, multiplying by |V p\ q ~~ 2 V p and 
then integrating over M 3 , we get, by Young's inequality, 



^IIVPII V, < 



\V p\ q \divv\dx + / p|Vp| 9_1 |Vdivv|dx 

+ 1 \Vv\\Vp\ q dx- - [ vV\Vp\ q dx 
Q Jm.3 



< llVplH, (||Vt;||ioc + ||p|| L oo||Vdivu||£«; 

H — / divw|V p\ q dx + ||p||z,°° ||Vdivv|| 
Q Jm.3 

< C\\V p\\ q Lq \\v\\ W 2, q + \\p\\ L oo\\Vdivv\\Li, 



(3.3) 
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since Ty 1,,? (R 3 ) L°°(R 3 ) as q > 3. Using Gronwall's inequality, we conclude that 



(3.4) 



l|Vp(t)||^ {R3 ) < ^||Vpo|||, + J ||p||L°°||Vdivu|| L gds^ exp (^J \\v\\ W 2, q ds 

< (llVpoIlL + V^lbHw) exp (Vt|M|w) • 
The proof is complete. □ 



3.2. Solvability of the deformation gradient with a fixed velocity. Due to the 

hyperbolic structure of (|3. lcj) . we can apply Lemma 13.11 again to solve the deformation 
gradient F in terms of the given velocity. For this purpose, we have 

Lemma 3.3. Under the same conditions as Theorem there is a unique function 

F :=T{v) G W 1 ' 2 (0,T;L 9 (M 3 )nL 2 (M 3 ))nL oo (0,T;VF 1 ' 9 (M 3 )nT^ ll2 (]R 3 )) 

which satisfies the equation (|3. lc|) . Moreover, the deformation gradient satisfies 

ll F llL°°(o,T ; iy 1 >'j(R3))nH 1 (K 3 ) - ( H F (°) II w 1 ^ (K 3 )ni? 1 (K 3 ) + V^ll v II wj dXp(VT\\v\\w) ■ 

Proof. First, we assume that v G C 1 (0, T; Cq°(]R 3 )), F G C$°(R 3 )). Then, we can 
rewrite (|3. lcj) in the component form of columns as 

d t F k + v-VF k = VvF k , for all 1 < k < 3. 

Applying Lemma 1331 with Aj(x,t) = Vj(x,t)I for 1 < j < 3, B{x,t) = —Vv, and f(x,t) = 
0, we get a solution 

oo 

F G pi {C^T,^ 1 ^ 3 )) nC(0,T;H s {R 3 ))}. 

s=3 

This implies, by the Sobolev imbedding theorems, 

oo 

F G pi C^O^-C^R 3 )) = C^T^QR 3 )). 

k=l 

Next, for v G W, by an argument of dense sets, there is a sequence of functions v n in the 
space CHO^C^QR 3 )), v n u in W, and Fft G Cg°(M 3 ), Fft F in ^'"(M 3 ) n H\R 3 ). 
Hence, v n — > v in C(S(0,a) x (0,T)) for all a > where 5(0, a) denotes the ball with 
radius a and centered at the origin. According to the previous result, there are {F n }^ =1 
satisfying 

d t F n + v n ■ VF n = Vv n F n , (3.5) 

with F n (0) = Fq, F n G C 1 (0, T; C°°(IR 3 )). Multiply ([33]) by |F n |P" 2 F n for any p > 2, and 
integrating over M 3 , by integration by parts, we obtain, using Young's inequality, 



™ / |F„|'V.r 
p dt 



— [ v n -V\F n \ p dx+ [ Vv n \F n \ p ~ 2 F 2 n dx 

P 7R3 J R 3 



, 1 + p,,„ IIP MY7 II 

p 
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Then, by Gronwall's inequality, one has 

/ \?n\ p dx(t)< [ \F n (0)\ p dxexp ( [ (p+l)||Vu n || L ~d 5 
Jm.3 Js? \Jo / 

< J |F n (0)| p cfeexp (J (p + l)\\v n \\ W 2, P ds\ . 

Thus, 

ll F n||L°°(0,T;LP(R3)) - ex P f ~ IbnllwV^ ||F„(0) ||i,p(r3) 
<exp (2||« n || w Vt) ||F„(0)|| Lpm3) < oo. 

Letting p — > oo, one obtains 

ll F n|U«(Q T ) < exp(||u„|| w \/t)||F n (0)|| L oo (R 3) < expdl^llw^llFn^H^i.^KS) < oo. 

Hence, up to a subsequence, we can assume that v n were chosen so that 

F„^F weak-* in L°°(0, T; L q (R 3 )). 

Taking the gradient in both sides of (|3.5|) . multiplying by |VF n | 9_2 VF n and then inte- 
grating over M 3 , we get, with the help of Holder's inequality and Young's inequality, 

VF 9 

q dt L 



(3.6) 



<2 f \VF n \ q \Vv n \dx+ [ |F n ||VF n | 9 " 1 |VVu n |(i3; - - / v n V\VF n \ q dx 
<C \VF n \ q \Vv n \dx + |F n ||VF n | 9-1 |Wv n |da; 

< C||VF n ||| 9 ||u n ||^2, 9 + ||F n ||x,c» ||u n ||v^2,g|| VF^H^g 

< C||VF n ||| 9 ||u n || W 2, 9 + C||t; n || l y2,g||VF n ||^ 1 , 
since q > 3. Using Gronwall's inequality, we conclude that 

||VF n (t)|| w(M 3) < M|VF„(0)||i? + J \\v n \\ W 2, q ds\ exp(J \\v n \\ W 2, q ds) 
< (||VF„(0)|| w + V~t\\v n \\w) exp{Vi\\v n \\ w ), 

and hence, 

II^F|| LO o( 0)T;L g( R 3)) < lirninf || VF n 11^=0(0^^9) 

< (||VF(0)|| L , + v^Hullw) exp(v / T||v|| w ). 
Thus, 

II f ||lcx.(o i t;W 1 .8(K 3 )) < (ll F (0)||i^i l9 + v^lMlw) exp(\/T||u|| w ) < oo. 

Passing to the limit as n — > oo in (13. 5p . we show that (j3. lcj) holds at least in the sense 
of distributions. Therefore, d t F £ L 2 (0, T; L 2 (M 3 )), then F e W 1,2 (0, T; L q (R 3 ) n L 2 (M 3 )). 
The proof is complete. □ 



(3.7) 
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3.3. Local solvability of (|3.1b[) . For simplicity of the presentation, we consider the case 
(j, = 1 and A = without loss of generality. In order to solve (j3.1b|) . we consider the 
bilinear form E(u, ip) and linear functional L(ip) defined by 

E(u, i/))= f (pd t u - Au - Vdivu, d t ip - k(Aip + Vdivip))dt 
Jo 

- (u(0),AV>(0) +VdivV>(0)), 

i-t 

L(ip) = - l (pv ■ Vv + VP — div(pFF T ), d t ip - k(Aip + Vdwip))dt 
Jo 

- (u ,A^(0) +Vdiv^(0)), 

with 

k = (2||p||l°°(q t )) _1 

for ifi S VV', where (•, •) denotes the inner product in L 2 . 

We first notice that L(ip) is a linear continuous functional of ip with respect to the norm 
Moreover, we have 

(WVpdt^Wh + k W A ^ + Vdiv^||| 2 - k(pd t ip, Aip + Vdiv^) dt 

+ \ (l|V^(T)|l! 2 + ||W(0)||| 2 + ||divV(T)||! 2 + ||div^(0)|| 2 L2 ) 

f T / 3 2k \ 

J Q (llvWIli* +fc||AV + Vdiv^||| 2 - -\\^pd t i,f L2 - y||AV + Vdiv^||| 2 J dt 

+ \ (HWCnill. + l|VV(o)||| 2 + ||div^(r)||| 2 + ||divv(0)||i 2 ) 

> c 1 1 -0 1 1 VV' ' 
for some Co > 0, since 

||A^ + VdivV|| L 2 > c^W\ m 

from the theory of elliptic operators. Hence, by the Lax-Milgram theorem (see [5]), there 
exists a u G VV' such that 

E(u,iP)=L&) (3.8) 

for every ip G VV'. 

Now, let ifi be the solution of the problem 

d t i> - k(A^ + Vdiv^) = 0, 

?(0) = h(x), 

with h(x) smooth enough. Replacing in (|3.8p ip by ip, one obtain 

(u(0) - u , Ah + Vdiv/t) = 0, 
which implies u(0) = Uo- Next, let tp be a solution of the following problem 

dti> - k(A4> + Vdiv^) = g(x, t), 

^(0) = 0, 



> 
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with g smooth enough. Replacing ip by %p in (|3.8p . one obtain 



/ (p<9 4 u - Au 
J o 



Vdivu + pt; • Vv - div(pFF T ) + VP, = 0. 



This implies that (u, p, F) satisfies (|3.ip a.e. in (0,T) x M 3 . 

Next, we prove the higher regularity for u; that is, u £ W. First, we multiply (|3.1b|) 
by <9fU, and use integration by parts and Young's inequality to obtain 



2 
3 



jf ||Vp9 t u||| 2 d S + ^||Vu(i)||£ 2 + |||div U (t)||| 2 
< |||Vu(0)||| 2 + i||divu(0)||| 2 - j\vP,d t u)ds+ 



t 

+ C sup \\p\\l°° I \\v\\ 2 L 6 



se(o,t) Jo 



Vv\\ 2 L3 ds+ [ (div(pFF T ), dtu)ds + -j / ||S t Vu||£ 2 ds 
< |[|Vu(0)||| 2 + ~||divu(0)||£ 2 + j* \\Vp\\l 2 ds + \f^ \\d t VM\\ 2 L ,ds+ 



— ( ' sup / ||«||x6||Vv||^ 3 ds + J ||VF||| 2 cZs + - J ||v^^ u lli3^' s - 



sG(0,i) JO 



(3.9) 

In particular, if we multiply (|3.1b|) by c^u, integrate over R 3 , and let t = 0, we obtain 



11^(0)^(0)11^ = / Au(0) + Vdivu(O) - p(0)v(0) • Vv(0) - VP(p(0)) 

+ div(p(0)F(0)F T (0))) d t u(0)dx 
< ||7^bT||Zi(l|Au || L2 + ||p(0)||r- ||«(0)|U<» || Vv(0)\\f 

+ C||Vp(0)|| i2 + ||VF(0)|| i2 )|| v ^(0)^u(0)|| z2; 



which implies that 

||Vp(0)^u(0)|| l2 < C(||Au || L 2 + IKIIlcoHVuoII^ + CHVyOoll^ + l|VF || L2 ). (3.10) 
Now we differentiate (|3.1b|) with respect to t so that we get 
dtpdtu + pdftXi — Adtu — Vdiv(<9 4 u) 



-dtpv ■ Vf — pdtv ■ Vv — pv ■ Vdtv — Vc^P + <9idiv(pFF ). 



(3.11) 
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Multiplying (|3.1ip by dtu, integrating over M 3 , and using the continuity equation, we 
obtain 

--\\^pd t u\\ 2 L2 + - / d tP \d t u\ 2 dx + ||VcH|| 2 + ||diva iU ||2 2 

2 at 2 J R 3 

< ||/>||L«||v||£=o||V«||ia||Vu||x3||^u|| L 6 + IHIloo \\v\\ 2 L6 \\V(Vv)\\ L 2\\d t u\\ L e 
+ \\Vp\\L^\\^/pdtv\\ L 2\\Vv\\ La \\d t \i\\ L 6 + ||p|U-|| , y|U-lklU6||Vi;|| I/ 3||V5tu|| I ,2 
+ \\y/ph^\\Vd t v\\ L 2\\v\\ L oo\\^pd t u\\ L 2 + \\d t P\\ L 4Vd t u\\ L 2 ( - ! - 2) 

+ \\Vd t u\\ L2 {\\Vp\\ L 2\\v\\ L o \\F\\ 2 L oo + \\p\\l™\\Vv\\ L 2\\F\\ 2 loo 

+ \\p\\l°°\\f\\l°°\\v\\l°°\\vf\\ L 2 ). 



Integrating (]3. 12|) with respect to t, using the continuity equation and the Gagliardo- 
Nirenberg inequality 

||Vu|||a < C||VuM|Au|| i2j 
we find, since p G L°°((0,T) x M 3 ), 

l||V^u(i)||i 2 +^(||V9 t u||| 2 + ||div5 t u||| 2 )d S 

< l\\Vf>d t u(0)\\h 

+ £ \\Vp\\l°° (llv / P^IU 2 ll^||| 2 ||Au||| 2 ||V5 t u|| i2 + |M| ioo ||V3 t u|| i2 ||Vpd t u[| L 2)d s 

+ / < ||Vp||L-||V^|| i2 || v || L ^||^u|| i2( i s + / ||^P|| L2 ||V9 t u|| i2 d s 

+ Ci / ||V^u||ia(||«||ioo + ||Vw||i3)ds 
Jo 

<i||V^ t u(0)||! 2 +C 5 ^ IblllccHv^ull^cfa + Ci (|| W ||l c||V W || 3 2 ||At;|| L2 

+ C7 5 ||Vt;||£ 2 ||A W ||2 2 ) ( fe + C (5 /"lIv^twIliallVrlliallAuH^cfa + Cs / Hdiv^H^ds 

JO Jo 

+ 5 [\\\Vd t u\\ 2 L2 + \\^d t u\\l 2 )ds + C 6 [\\\v\\ 2 Lao + \\Vv\\ 2 L2 )ds+ [ \\Vd t v\\ 2 L2 ds, 
Jo Jo Jo 

(3.13) 

where 5 is a small constant, C\ depends on v, po, and Fo, since 

Vp£ L°°(0,T;L 9 (lR 3 )nL 2 (M 3 )), VF G L°°(0, T; L 9 (IR 3 ) n L 2 (R 3 )) 

as stated in Lemma [3. 2 1 and Lemma [3 .31 Since we are only interested in the local existence, 
so we can restrict t < T < 1. 
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Adding (|3.9|) and (|3. 13|) for some suitable S < i, one obtains, first by Gronwall's 
inequality, 

ll\/P^t u lli oo (0,t;i 2 (R 3 )) - C; 

and, second, 

\\\VpdMt)\\h + I j\\\Vd t u\\ 2 L2 + ||diva t u||| 2 )d S + ~ jf ||^u||| 2 d S 

+ ^l|Vu(t)||| 2 + i||divu(t)||| 2 <C7, 
which implies 

f ^ ugL °°(°' t ; l2 ( m3 )); ftue^o.T;]? 1 

I Vu G L°°(0,T;L : 

On the other hand, we can rewrite (|3.1bp as 

-Au - Vdivu = -pd t u - pv ■ Vv - VP + div(pFF T ), 

which is a strongly elliptic equation. Hence, by the classical theory in the elliptic system, 
we deduce that 

u G L°°(0,T;# 2 (R 3 )), 
since from Lemmas I3.2H3.3I and ()3.14j) one has 

-pd t u — pv • Vv — VP + div(/)FF T ) G L°°(0, T; L 2 (M 3 )). 

Moreover, since p is bounded from below and y/pd t u G L 2 (0, T; L 2 (M 3 )), we know that 
9tU G L 2 (0, T; L 2 (M 3 )). Hence, by the Gagliardo-Nirenberg inequality, as q G (3,6], 

ll^t u llL 2 (Q,T;L9(R3)) < C||^ u IIl 2 ((O,T)xR 3 )II V ^ u IIl2( ( O>T ) x1 ;3), 

for some 9 G [0,1). This implies that, by (|3TTi|) . <9iU G L 2 (0, T; L 9 (M 3 )). Thus, by the 
classical elliptic theory, we obtain u G L 2 (0, T; W 2 '«(R 3 )) since now 

-pd t u - pv • Vv - VP + div(pFF T ) G L 2 (0, T; L«(M 3 )). 

Hence, we can conclude that u G W. 

3.4. Existence for The above argument leads us to define the map 

u = H(v) 

from W to itself through the maps g : v h-> 5(f), / : u h-> T(i>) and d : (5(t>), t>, T(v)) \-t u. 
Hence the solution of (jl.ip is obtained from a fixed point of the map T~L. To find a fixed 
point of we will use the Schauder-Tychonoff fixed point theorem (Theorem 5.28, |28]). 
Define 

M = {ip : max(||V'||x2(o,T;W2,<r(R3) niJ -2( R 3)), || VpdtV'lU^o.T^QRS)), 

ll^lll,«=(0,r;i?2(R3)), ||3^||l2(0,T;H1(I: 3 ))) ^ r }' 

where 

r = C (||u Q || H 2 + HuollioollVuollia + ||^ollvF 1 .9nH 1 + H F ollvK 1 .9nH 1 ) 
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with some sufficiently large Cq > 0. Clearly, M is a compact and convex set in L 2 ((0, T) x 
M 3 ). Hence, we need to show that U{M) C M (i.e., U maps M into M) and % is 
continuous in M with respect to the norm in L 2 ((0,T) x R 3 ). 

We first prove that T-L(M) C M for some T = T. Indeed, assuming v £ M, from Lemma 
T2l and Lemma 13.31 we know that 



aexp(— ry/i) < p(x,t) < /3exp(ry/t); 

IpIU^o^W^qr 3 )) < (IIpoIIwl^R 3 ) + Vtr) exp (Vtr) ; (3.15) 
l F IU<»(0AW rl .«CR 8 )) - (ll F ( )llw 1 .'J(M^) + V^r) exp( v / tr). 
Hence, from (|3.13p and (|3.15p . it follows that 



l\\VpdMt)\\l 2 + ~J*\\Vd t uf L2 ds 



< \\\^pd t u(0)\\ 2 L2 + C 5 t(r 2 + ^[IV^uHloo^^s)) + C s t(r A + r 6 ) + Ctr 2 . 

Using (|3.10p . and taking 5 and T sufficiently small, we derive from the above inequality 
that 

II\Z^ U IIl°°(0,T;Z,2(k3)) + H V ^ U Hl 2 (0,T;L 2 (K3)) - ^ ' 

Since p is bounded from below, we obtain 

2 2 ^ 2 

IIV^ U IIl°°(0,T;L 2 (]R 3 )) + H^ tU llL2(0,T;H 1 (R 3 )) - Q F ' (3.16) 



To estimate the norm HuH^oomTj^fR 3 ))) noticing that 

-Au - Vdivu = -pd t u - pv-Vv- VP(p) + div(/?FF T ), 
one has from the theory for elliptic equations that 

II U IIl°°(0,T;H 2 (R 3 )) - llv^lli 00 II\/P^ u IIl°°(0,T;L 2 ) + \\p\\l°° II ^ v IIl°°(0,T;L 2 ) 

+ cIIVpIIl^o^l 2 ) + II vf IIl°o(o,t ; l 2 ); 

and hence 

||Au + Vdivu[|£ M(0j5 ;. H2) < r 2 . 

Next, we need to obtain estimates on ||u.||^2 (q ) t , ;PF 2 '*(K 3 )) ■ Indeed, we have, by the 
classical theory of elliptic equations, 

/ HuH^dt < C / (HvHioo ||V«||i, + ||Vp||i, + \\d t u\\ 2 Lq + || VF[|£,)di 
J io 

< Cr A T + r 2 T< r 2 , 

for some sufficiently small T. Hence, we show that T~L{M) C M. 

Finally, we need to prove the continuity of T~L in M. First we observe that if {v n }^ =1 C 
M, then there exists a subsequence (still denoted by {v n }^ =1 ) such that as n — > 00, v n — > v 
strongly in M . Let p n and p be the solutions of 

d t pn + div(p n v n ) = 0, 

and 

d t p + div(pv) = 0, 
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with p n (0) = p(0) = po, respectively. Denoting p n = p n — p, then p n satisfies 

dtPn~ + v n ■ Vpn~ + (v n -v)-Vp + p^divv n + pdiv(v n - v) = 0, 
with /5n(0) = 0. Repeating the argument in Lemma 13.21 we have 

||^|| 2 2 < exp(CrT) / (\\(v - v n ) ■ V p\\\ 2 + \\pd\w{v n - v)\\ 2 L2 ) dt, 
Jo 

which implies that p n — > p strongly in L°°(0, T; L 2 (1R 3 )). Similarly, we can show that F n — > 
F strongly in L°°(0, T; L 2 (M 3 )). Now let u ra and u be the solutions of (]3.1b[) corresponding 
to v n and v with u n (0) = u(0) = uo respectively. Then one has, denoting U n = u n — u 
and V n = v n — v, 



Multiplying the above equation by dtU n , integrating over (0, T) x M 3 , and thanking to 
the convergence of p n and F n , we can prove as a routine matter that VU n — > strongly 
in L 2 ((0,T) x M 3 ) and ^pn~d t U n strongly in L 2 ((0,T) x M 3 ). Due to the lower 
bound of p n , we deduce that dtU n — > strongly in L 2 ((0, T) x R 3 ), and hence U n — > in 
L 2 (0,T) x K 3 ) by using the identiy U n (t) = f*d t U n ds since U n (0) = 0. Thus, the map U 
is continuous in M. The existence of a local solution is completely proved. 



In this section, we will prove the uniqueness of the solution obtained in the previous 
section. Notice that, the argument in Section 3 yields that 



d t u E L 2 (0, T; L 2 n L q (M. 3 ), Vp G L 2 (0, T; L 2 n L 9 (M 3 )), VF € L 2 (0, T; L 2 n L q (R 3 )), 



Pnd t U n - AC/ n - VdivU n 



+ div(p n F n F r l - p¥¥ 1 ) - VP(p n ) + VP(p). 



4. Uniqueness 



for g > 3. Hence, using the interpolation, we see that 




(4.1) 
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which yields 
d_ 
~dt 



Hilars) < ||divu 1 || L ^||r||| 2 +e||Vt;||| 2 +C7(e)||V5(u 2 )r||2 f 



+ e\\Vv\\i 
< ||divui|| Ll 1 ••' 1 



i2 +C( £ )||5(u 2 )||| c 



l 2 + e\\Vv\\i 2 + C(e)\\VS{u 2 )\\i z \\r\\l 2 



WW 2 

\ r \\2 



+ e||Vt;||^ + C( e )||5(u 2 )||i 

where m (e) = ||div Ul || L ~ + C(e)(\\ VS(u 2 )||£ s + ||S(u 2 )||Lo). 
Similarly, from (|3.1c|) . we obtain 

f d t G + ui • VG + v ■ VT(u 2 ) = VuiG + WT(u 2 ), 
( G(0) = 0. 



Multiplying 



3]) by G, and integrating over 



, we get 



which yields 
It 



\\G\\ 



|G| 2 divuidx + [ v- VT(u 2 ) : Gdx 
G T V Ul G(ij; + / VvT(u 2 ) : Gdx, 

< ||divui|| ioo ||G||| 2 +e\\Vvf L2 + C(e)||VT(u 2 )G|| 



+ e||Vw||i,+C(e)||r(u a )||i-|^||l 
< ||divu 1 ||^||G||| 2 +e||V W ||| 2 +C7( e )||Vr(u 2 )||| 3 ||G||| 2 

+ £ ||VHH 2 +C(e)||r(u 2 )||| 0O ||G||i 
<r ?2 ( e )||G||| 2 +2 £ ||Vi;||! 2 , 

where 773(e) = ||divui||i- + C(e)(||VT(u^ + \\T(u 2 )\\ 2 Lx ). 
For each Uj j = 1,2, we deduce from (|3.1bj) that 

S(uj)dtUj — fiAuj — (fi + A)Vdivuj 

= -5( Uj )( Uj • V) Ui - VP(5( Uj )) + div(5(u J )r(u J )T(u J ) T ), 
Uj(0) = u . 

Subtracting these equations, we obtain, 

S(ui)d t ui — S(v 2 )dtU2 — /uAw — (p + A)Vdivt> 
= -5(ui)(ui • V)ui + cS(u 2 )(u 2 • V)u 2 - VP(5(ui)) + VP(«S(u 2 )) 
+ div(5(u 1 )T(u 1 )T(u 1 ) T ) - div(S(u 2 )T(u 2 )T(u 2 ) T ). 



Since 



<S(ui)(ui • V)ui + <S(u 2 )(u 2 • V)u 2 

: -S(ui)(v ■ V) Ul - (5(ui) - S(u 2 ))(u 2 • V)ui - 5(u 2 )(u 2 • V)v, 
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and 

<S(u 1 )T(u 1 )r(u 1 ) T - S(u 2 )T(u 2 )T(u 2 ) T 
= S(u l )GT(u 1 )) T + rr(u 2 )r( Ul ) T + S(u 2 )T(u 2 )G T , 
we can rewrite (14, 5ft as 
S(ui)dtv — /jAv — (/i + A)Vdivv 

= -rd t u 2 - 5(ui)(« • V)ui - r(u 2 • V)uj - 5(u 2 )(u 2 • V)v - VP(5(m)) + VP(«S(u 2 )) 
+ div(cS(u 1 )GT(u 1 ) T + rT(u 2 )T( Ul ) T + S(u 2 )T(u 2 )G T ). 

(4.6) 

Multiplying (j4.6j) by using the continuity equation (|l.la|) and integrating over M 3 , we 
deduce that 

~ / S(u 1 )\v\ 2 dx+ [ (p\Vv\ 2 + (X + fi)\divv\ 2 )dx 

= \ -5(ui)(ui • V)v ■ v — rdtU2V — S(u\)(v ■ V)u\v — r(u 2 • V)uiu 

- 5(u 2 )(u 2 • V)vv - VP(5(uO)« + VP(S(u 2 ))v 

- (5( Ul )GT( Ul ) T + rT(u 2 )T(u 1 ) T +5(u 2 )r(u 2 )G T )V7;} ( ix 

< ellVvllia +C(e)[|5(u 1 )||i« a o|[ixi]|i«»[|«[|i a +e]|V«|[i a +C(e)||a t ii2[|is[k[l! #a (4.7) 
+ ||5(ui)|Uoo[|Vui|| x <»||t;||i a +2||u 2 ||L«»||Vui||iao([|r||i a + \\v\\ 2 L2 ) 
+ e\\ Vvf L2 + C(e)[|5(u 2 )||ioc \\u 2 \\ 2 Lao |M|| 2 + e\\Vvf L2 
+ C( £ )(sup{P'(7 ? ) : C{T)^ < r, < G(T)}) 2 ||r||| 2 + e[|Vw[|£ a 
+ C( £ )(|| 4 S(u 1 )||| 00 ||r(u 1 )||| 00 ||G|| 2 L2 

+ ||5(u 2 )|| 2 LOO [ircuajiii-iiGii^ + [irii^iircuOHLooiircua)^. 

<5e||V.;|| 2 2 +7 ? 3(e)(||r||| 2 + H 2 , 2 + ||G|| 2 2 ) 

with 

773(e) = C(e)\\S(u 1 )\\ 2 Loo \\u 1 \\l ao +C(e)\\d t u 2 \\ 2 L:i + ||5(ui)|U«, [|Vui|| £ oo 
+ 2||u 2 ||i=c||Vu 1 ||ioo +C(e)||5(u2)||ioo||u 2 ||ioo 
+ C(e)(sup{P'(r / ) : G(T) _1 < r, < C(T)}) 2 

+ C( e )([|5(ui)[|i B o[|r(ui)||i. + [|5(u a )||i«, ||T(u 2 )[|i.) 
+ l|r(ui)[|i fl o||T(u a )[|i B o. 
Summing up (|4.2p . (|4.4p . and (|4.7p . by taking e = ^, we obtain 

4/ (5^i)|«| 2 + |r| 2 + |G| 2 )dx + / L i / IV^I 2 ^ 
«t 7r3 Jrs 

< 2(77 3 ( £ ) + 772(e) + r ?1 (e))(|| V ||i 2 + ||r||| 2 + ||G|| 2 2 ) (4.8) 

<2r/(e,t) / (S( Ul )H 2 + |r| 2 + |G| 2 )dx, 
Jr 3 
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with 



min{min x6R 3 S(u 1 )(x,t), 1} ' 
It is a routine matter to establish the integr ability with respect to t of the function rj(e,t) on 
the interval (0, T). This is a consequence of the regularity of ui, 112 € W and the estimates 
in Lemmas 13.21 and 13,31 for 7~(uj) with i = 1, 2. Therefore, ()4.8|) . combining with 

Gronwall's inequality, implies 



/ 



(5(ui)|v| 2 + |r| 2 + \G\ 2 )dx = 0, for all te(0,T), 
and consequently 

u = 0, r = 0, G = 0. 
Thus, the proof of uniqueness is complete. 
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